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Abstract 



The temperature dependence of the mass, leptonic decay constant, and width of 
heavy-light quark peseudoscalar and vector mesons is obtained in the framework 
of thermal Hilbert moment QCD sum rules. The leptonic decay constant of both 
pseudoscalar and vector mesons decreases with increasing T, and vanishes at a 
critical temperature T c , while the mesons develop a width which increases dra- 
matically and diverges at T c , where T c is the temperature for chiral-symmetry 
restoration/quark-gluon deconfinement. These results indicate the disappear- 
ance of hadrons from the spectral function, which then becomes a smooth func- 
tion of the energy. This is interpreted as a signal for deconfinement at T = T c . 
In contrast, the masses show little dependence on the temperature, except very 
close to T c , where the pseudoscalar meson masses increase slightly by 10-20 %, 
and the vector meson masses decrease by some 20-30 %. 
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1 Introduction 



The thermal behaviour of hadronic Green's functions, obtainable in a variety of 
theoretical frameworks, plays a fundamental role in understanding the dynamics 
of the quark-gluon plasma. One such framework is that of QCD sum rules pQ, 
based on the Operator Product Expansion (OPE) of current correlators beyond 
perturbation theory, and on the notion of quark-hadron duality. The extension 
of this program to finite temperature was first discussed long ago in pj. It is 
based on two basic assumptions, (a) that the OPE remains valid at T ^ 0, 
with perturbative QCD (PQCD) and the vacuum condensates developing a 
temperature dependence, and (b) that the notion of quark-hadron duality also 
remains valid. Additional evidence, from solvable quantum field theory models, 
supporting these assumptions was provided later in [3] . Numerous applications 
of this technique have been made over the years [5], leading to the following 
consistent picture, (i) With increasing temperature, particles that are stable at 
T = develop a non-zero width, and resonances become broader, diverging at 
a critical temperature interpreted as the deconfinement temperature (T c ). This 
width is a result of particle absorption in the thermal bath, (ii) The onset of the 
continuum in hadronic spectral functions, traditionally accounted for by PQCD, 
decreases and approaches threshold near T c . In other words, as T — > T c hadrons 
melt and disappear from the hadronic spectral functions, which become perfectly 
smooth, (iii) This scenario is further supported by results for hadronic and 
electromagnetic mean-squared radii, which also diverge at T c [3]. In addition, 
QCD sum rules in the axial-vector channel have provided (analytical) evidence 
for the near equality of the critical temperatures for deconfinement and chiral- 
symmetry restoration [5] . Regarding the temperature dependence of a hadronic 
mass, this parameter does not appear to be a relevant signal for deconfinement. 
Conceptually, given either the emergence or the broadening of an existing width, 
together with its divergence at T c , the concept of mass looses its meaning. In 
practical applications, in some cases the mass increases slightly with increasing 
T, and in others it decreases. 

In this paper we use Hilbcrt moment QCD sum rules for heavy-light quark 
pseudoscalar and vector meson correlators to determine the thermal behaviour 
of the hadronic masses, couplings, and widths. At T — this problem was 
discussed long ago in [7j-[S]. While there are only four ground-state pseudoscalar 
heavy-light quark mesons in the spectrum (£>, D s , B, and B s ), and similarly 
for vector mesons, it is possible to determine the decay constants for arbitrary 
meson masses in a self-consistent way [5]. The result is that the leptonic decay 
constants obey a scaling law as a function of the meson mass. Since we are only 
interested in the temperature behaviour of the hadronic parameters, we shall 
normalize all our results to those at T = 0, thus obtaining a universal functional 
relation. We find that the meson masses are basically independent of T, except 
very close to T c where they increase slightly (pseudoscalars) by 10 - 20 %, or 
decrease (vector mesons) by 20-30 %. Here T c is the critical temperature for 
chiral-symmetry restoration/quark-gluon deconfinement. The leptonic decay 
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constants decrease with increasing T, and vanish at the critical temperature. 
Pscudoscalar mesons develop a non-zero hadronic width that increases with 
increasing T and diverges at T c , while vector meson widths exhibit a similar 
beahaviour. These results may be interpreted as providing (analytical) evidence 
for quark deconfinement at T = T c . 



2 Pseudoscalar mesons 



We begin by defining the correlator of axial- vector divergences at finite temper- 
ature, i.e. the retarded Green's function 



ip 5 (q 2 ,T)=i / d 4 xe^ x 6(x ) « \ [WA^x) , d v A\ 



» , 



(1) 



where d^A^x) = ttiq : q(x) Q(x) : , q (Q) refers to the light (heavy) quark 
field, and vtlq » m q is assumed. The matrix element above is the Gibbs 
average 



« A ■ B »- 



exp(-E n /T) < n\A ■ B\n > /Tr{exp{-H/T)) 



(2) 



where \n > is any complete set of eigenstates of the (QCD) Hamiltonian. We 
use here the quark-gluon basis, which allows for a smooth extension of the QCD 
sum rule program to non-zero temperature 3 . At T = and to leading order 
in PQCD 



• Imip 5 (x, 0) 



{1-xf 



(3) 



where x = itiq/s, with s > rriq, and < x < 1. At finite temperature there are 
two distinct contributions to the correlator, the so called scattering term (q 2 
space- like), and the annihilation term (q 2 time- like) [2]. After a straightforward 
calculation we find the former to be exponentially suppressed, so that it can be 
safely neglected, while the latter is given by 



1 Imip 5 (x,T) = \hli ( 'r, !.r. (I j <jl-//p — (i -f.r) 



■Up 



2T 



0--X) , (4) 
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where Imip5(x,0) is given by Eq.(3), np(z) = (1 + e 2 ) -1 is the Fermi thermal 
function, and in the rest frame (q = 0) of the thermal bath x = mq/uj 2 . 
At temperatures below T ~ 200 MeV one can safely assume the heavy quark 
mass to be temperature independent [TUj. The first thermal function above is 
exponentially suppressed and can be safely neglected for temperatures of order 
0(100 — 200 MeV) , but the second one does contribute near threshold. 
Up to dimension d — 6 the non-perturbative expansion of the correlator at 
T = is given by [7], [3] 



fp5(q 2 )\NP 



m Q m Q T 

2 7T C*4 < O4 > H — ~. 7; TT^T C5 < O5 > 

4 (rrip — q^) 6 



m" Q -q 
2 



ii>7 



111; 



{m 2 Q - q 2 ) 2 (m 2 Q - q 2 ) 3 {m 2 Q - q 2 ) 4 
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C 6 < 6 >, (5) 



where 



C 4 < O4 >= — — < a s G 2 > - m Q < qq>, 



(G) 



C 5 < O5 >=< g s qi<T^ G a \ a q>=2M 2 <qq> , 



(7) 



C 6 < O e >= 



7ra, 



vs 16 . l2 /n x 

q> =4> -—-Ka s p\ <qq> r , (8) 



and m c (m c ) ~ 1.3 GeV, to 6 (to 6 ) ~ 4.2 GeV, < ^ >~ (—250 MeV) 3 , 
< a s G 2 /12vr >~ 0.003 - 0.006 GeV 4 , Af 2 ~ 0.4 - 0.6 GeV 2 , and p ~ 3 - 5 
accounts for deviations from vacuum saturation (VS). Use of these values in 
Hubert moment sum rules reproduce the pseudoscalar meson masses at T = 0. 
However, they will play no crucial role in our analysis as changes in these pa- 
rameters would only affect the normalization at T = 0. 



For the light-quark condensate at finite temperature we shall use the result of 
|llj . obtained in the framework of the composite operator formalism, valid for 
the whole range of temperatures: T = — T Cl where T c is the critical tempera- 
ture for chiral symmetry restoration. There is lattice QCD evidence [12] as well 
as analytical evidence [6] for this critical temperature to be (almost) the same 
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as that for deconfinement. The ratio R(T) =« qq » / < qq > from [IT] as 
a function of T/T c is shown in Fig. 1 . 



The low temperature expansion of the gluon condensate is proportional to the 
trace of the energy-momentum tensor, and it starts only at order T 8 |13j . viz. 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

T/Tc 



Figure 1: The light-quark condensate ratio R{T) =<< qq >> / < qq > as a 
function of T/T c from p]. 

Because of this T- dependence, the gluon condensate remains essentially con- 
stant up to temperatures very close to T c . Hence, the thermal non-perturbative 
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QCD correlator is basically driven by the quark condensate. Concerning the 
dimension d = 6 condensate, it has been argued that the vacuum saturation 
approximation breaks down at finite temperature [14] . This is based on the 
comparison between the slopes of the low temperature expansion (chiral pertur- 
bation theory) with and without assuming vacuum saturation. They are in fact 
numerically different. However, this result is only valid at very low temperatures 
(T << /„.); hence it cannot be extrapolated to T ~ T c . In fact, both the quark 
condensate and the four-quark condensate should vanish at the same tempera- 
ture T = T c . In any case, numerically, at temperatures of order T ~ 100 Me V 
the quark condensate dominates over the gluon condensate, the dimension d = 5 
condensate is comparable to << qq >>, and the dimension d = 6 condensate is 
almost two orders of magnitude smaller. Hence, potential violations of vacuum 
saturation can be safely ignored. Finally, at finite temperature it is possible, 
in principle, to have non-zero values of non-diagonal (Lorentz non-invariant) 
vacuum condensates. There is one example discussed in the literature [15] with 
enough detail to make a numerical estimate of their importance, and it refers 
to operators of spin-two (quark and gluon energy momentum tensors) . The low 
temperature expansion of these terms starts at order 0(T 4 ), in contrast to a 
T 2 dependence for the diagonal condensates. We find that at temperatures of 
order T ~ 100 MeV both non-diagonal condensates are three orders of magni- 
tude smaller than the corresponding diagonal equivalents. We shall then ignore 
non-diagonal condensates in the sequel. 

Turning now to the hadronic sector, the spectral function at T = can be 
written as 



where Mp and fp are the mass and leptonic decay constant of the pseudoscalar 
meson, and the continuum, starting at some threshold so, is modeled by per- 
turbative QCD. With this normalization, f„ ~ 93 MeV. Anticipating the pseu- 
doscalar mesons to develop a sizable width Tp(T) at finite temperature (particle 
absorption in the thermal bath) , and using a Breit-Wigner parametrization, the 
following replacement will be understood 



where the mass and width are T— dependent, and the constant is fixed by re- 
quiring equality of areas, e.g. if the integration is in the interval (0 — oo) then 
const = 2MpT 'p/w. The continuum threshold sq above also depends on temper- 



- Imip 5 (s)\ H AD = 2fpM P S{s-Mp) + 0{s- So )-Imyj 5 {s)\p QC D , (11) 
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0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

T/T c 



Figure 2: The ratio So(T)/sq(0), Eq.(14), as a function of T/T c for tjiq = m c . 



ature; to a very good approximation it scales universally as the quark condensate 
HE], i.e. 



s (T) « qq » 
s (0) ~ < qq > 



where so(0) is clearly channel dependent. At the critical temperature we expect 
sq(T c ) = rriQ, in which case Eq. (13) can be rewritten as 



s (T) = « qq » 
S (0) ^ < qq > 

This is shown in Fig. 2 for the case toq = m c and sq(0) — 5GeV 2 ; a qualita- 
tively similar behaviour is obtained for itiq = and sq(0) ~ (1.1 — 1.3)Mg. 



«o(0) 



«o(0) 



(14) 
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T/T c 

Figure 3: The ratio M P (T)/M P (0) as a function of T/T c . 



The correlation function ip 5 (q 2 ,T), Eq.(l), satisfies a twice subtracted dispersion 
relation. To eliminate the subtractions one can use Hilbert moments at Q 2 = 
-q 2 = 0, i.e. 



(15) 

where N = 1, 2, .... Invoking quark- hadron duality 

<P {N \T)\ HAD = ^ n \T)\ QC d , (16) 



and combining the continuum contribution in the hadronic spectral function 
with the PQCD piece of the QCD counterpart leads to the finite energy Hilbert 
moments 
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0.9 1.0 



Figure 4: The ratio / P (T)// P (0) as a function of T/T c . 



i r o(T) i r o(T) rfs 

7mV'5(s,T)| P oL P = - Imip 5 (s,T)\ P QCD 

n Jm 2 Q s 

+ <pW{T)\ NP , (17) 



*N+2 



where Im ip5(s, T)\pole is given by the first term in Eq.(ll) modified according 
to Eq.(12), the PQCD spectral function corresponds to Eq.(4), and 



<pW{T)\ 



NP — —- 



m Q << qq » 



m 



2N+2 

Q 



1 _J_«q.Q»» _l (jV + 2)(jV + 1) 
YI-kuiq « qq » 4 V 



M 4 

— f - — (N + 2)(N 2 + 10 N + 9) 7T a s p 



« qq » 



rn 



Q 



81 



rn 



Q 



(18) 
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0.9 1.0 



Figure 5: The width T P (T)) as a function of T/T c , with r P (0) = 0. 



Using the first three moments one obtains the temperature dependence of the 
mass, the leptonic decay constant, and the width. Results from this procedure 
are shown in Figs. 3-5 for the charm case; in the case of beauty mesons, results 
are qualitatively similar. 



3 Vector mesons 

We consider the correlator of the heavy-light quark vector current 



U^(q 2 ,T) = i J d 4 xe^ x 9(x ) « \[V„(x) , Vj(0)]| » 

= -(g^q 2 -q^)H^(q 2 ,T) + q^ \q 2 ,T) , (19) 
where V^{x) =: q(xYf l _ l Q{x) :. In the sum rule analysis we shall use the function 
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— Q 2 H^(Q 2 ,T), which is free of kinematical singularities. A straightforward 
calculation gives 



-ImU^(x,T) = -i-(l-x) 2 (2 + x) 

7T 077 



1 - n F (z+) - n F (z-) 



(20) 



where z± = (1 ± x). In the hadronic sector, we define the vector meson 
leptonic decay constant f v through 



<0\Vp(0)\V(k) >=V2M v f v e tl , 



(21) 



so that the pole contribution to the hadronic spectral function is 2fyd(s — 
My). At T = the vector meson D*(2010) has a very small width in the keV 
range (96 ± 22 keV), which we expect to increase with increasing T, so that the 
replacement in Eq.(12) will be made. 




0.9 1.0 



Figure 6: The ratio M v (T)/M v (0) as a function of T/T c 
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Figure 7: The width T V (T) as a function of T/T c . 



The Hilbert moments at Q 2 = of the function -Q 2 nW (Q 2 , T) are given by 



(-) 



JV+l 



(JV + i)! VdQ 2 



ds 



) N+1 [-Q 2 n«(Q 2 ,r)]| Q2=0 



(22) 



Following the same procedure as for the pseudoscalar mesons (see Eq. (17)), 
the sum rules become 



1 r s o( T ~) Jo i /•sol-') w s 

i/ o ^/mII«( S ,T)| POiB = JfLlmnW(s,T)\ PQCD 



2 S 1 

Q 



(23) 
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where tp( N > (T)\np is given by 



NP 



mQ « qq » 



n 2N+4 



« a s G 2 » 
YI-kvtlq << qq >> 



(N + 2){N + 3) 



M 2 

x _°- + -l.(jV + 2)(20 +N -N' 2 )ira s p 

m Q 81 



<< qq » 



(24) 



Using the first three Hilbert moments to find the temperature dependence of 
the hadronic parameters, we obtain for the mass and the width of D*(2010) 
the results shown in Figs. 6-7. The behaviour of the vector-meson leptonic de- 
cay constant will not be shown, as it is essentially the same as that of the 
pseudoscalar-meson in Fig. 4. Similar results are found for the case of the beauty 
vector meson B*. 



4 Conclusions 



The thermal behaviour of pseudoscalar and vector meson decay constants, masses, 
and widths was obtained in the framework of Hilbert moment finite energy QCD 
sum rules. This behaviour is basically determined by the thermal light quark 
condensate on the QCD sector, and by the T-dependent continuum threshold 
on the hadronic sector. Normalizing to values at T = 0, and using the method 
of [5] for arbitrary masses, there follows a universal relation for the hadronic pa- 
rameters as a function of T/T c . Results show that the decay constants decrease 
with increasing temperature, vanishing at T — T c , while the widths increase 
and diverge at the critical temperature. Such a behaviour provides (analytical) 
evidence for quark-gluon deconfincment, and is in qualitative agreement with 
corresponding results obtained in the light-quark sector. Finally, pseudoscalar 
meson masses increase slightly with temperature by some 10 - 20 %, while vec- 
tor meson masses decrease by 20 - 30 %. Given the dramatic emergence of 
monotonically increasing widths r(T), there is little if any significance of this 
temperature behaviour of the masses, i.e. the relevant signals for deconfine- 
ment are the vanishing of the leptonic decay constants and the divergence of 
the widths at T = T c . 
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